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(N ; ABSTRACT 

^ , We report long-baseline interferometric measurements of circumstellar dust 

I i around massive evolved stars with the MIDI instrument on the Very Large Tele- 

CN I scope Interferometer and provide spectrally dispersed visibilities in the 8-13 fim 

I— li wavelength band. We also present diffraction-limited observations at 10.7 /xm on 

^ the Keck Telescope with baselines up to 8.7 m which explore larger scale struc- 
ture. We have resolved the dust shells around the late type WC stars WR 106 
and WR 95, and the enigmatic NaStl (formerly WR 122), suspected to have 

d \ recently evolved from a Luminous Blue Variable (LBV) stage. For AG Car, the 

I protoypical LBV in our sample, we marginally resolve structure close to the star, 

^ I distinct from the well-studied detached nebula. The dust shells around the two 

I WC stars show fairly constant size in the 8-13 fim MIDI band, with gaussian 

QQ I half-widths of ~ 25 to 40 mas, and the Keck observations reveal an additional 

CN I extended structure around WR 106. The visibility profiles for NaStl obtained 

^ I from two MIDI baselines indicate a compact source embedded in an extended 

I structure. The compact dust we detect around NaStl and AG Car favors recent 

^ I or ongoing dust formation. 

I Using the measured visibilities, we build spherically symmetric radiative 

^ ' transfer models of the WC dust shells which enable detailed comparison with 
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existing SED-based models. Our results indicate that the inner radii of the shells 
are within a few tens of AU from the stars. In addition, our models favor grain 
size distributions with large (~ 1 /xm) dust grains. This proximity of the inner 
dust to the hot central star emphasises the difficulty faced by current theories in 
forming dust in the hostile environment around WR stars. Although we detect 
no direct evidence for binarity for these objects, dust production in a coUiding- 
wind interface in a binary system is a feasible mechanism in WR systems under 
these conditions. 

Subject headings: Interferometry: general — Wolf-Rayet stars: individual(WR 
95, WR 106, NaStl) — Luminous Blue Variables: individual(AG Car) 



Introduction 



The most massive stars (> 4OM0) will likely pass through a Luminous Blue Variable 
(LBV) and Wolf-Rayet (WR) phase before ending their lives in a supernova explosion. LBVs 
and some WRs are known to show infrared excess associated with circumstellar dust. The 
LBV phase is shor t-lived and unstable with irre gular mass ejections of 10~^ to 1O~^M0 
for each outburst (IHumphreys fc Davidson! Il994l ) , resulting in circumstellar nebulae and 
associated dust. The dust in this phase is cons istent with grains of differing compo sitions 
ranging over metallic, silicate and carbon-rich (jVoors et al.l Il999l : IClark et al.ll2003l ). The 
WR stage, chara cterised by m assive winds (of the order M = lO~^M0yr~^) is divided into 
two main phases (ISmithlll968l ). the WN phase with strong spectral lines of N and CNO-cycle 
products, and the later WC phase with strong lines of C and other He-burning products. 
Dust formation is associated with the WC phase, consistent with the over-abundance of C, 
and is notabl y absent in the WN p hase (Ivan der Hucht et al.ll200ll ). although recent Spitzer 
observations (IBarniske et al.ll2006l ) may be the fi rst indications to the contrary. Long term 
and detailed infrared photometry of WC stars (jWilliams et al.l 119871 . hereafter WvdHT), 
has led to considerable progress in understanding the nature of the dust and its formation. 
The favored mechanism fo r dust-forrn ation in this stage is through wind-wind collision in 
a WR-OB binary system (lUsovl Il99ll ). F or the episodic dust produ cers WR 140 and WR 
137, photometric and spectrosocopic data ( Williams et al. 199ol . [200ll ) had already indicated 
bina rity which has been confirmed by long-baseline interferometric observations for WR 
140 (IMonnier et al.ll2004al ). High resolution images in the near-IR using aperture- masking 
techniques with the Keck telescope provide direct evide nce for binarity in the case of the 



const antly dusty, late-type WC stars WR 98a, WR 104 (IMonnier et al. 



19991 ). and more recently, the "cocoon stars" in the Quintuplet cluster (ITuthill et al 



I999I: iTuthi 



et al. 



20061). 
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Monnier et al.l (120071 ) present a strong case for binarity for all the dusty WRs in the Keck 



aperture-masking sample. However, whether all dusty late type WCs are b inaries is an open 



quest ion and dust formation in single WR stars has not been ruled out (jCherchneff et al. 



20001). 



In this paper, we present and interpret interferometric mid infrared observations of the 
late type WC stars WR 95 and WR 106, the suspected LBV-WR transition star NaStl, 
and the LBV AG Car. The main goals of this campaign were to resolve the sizes of the 
dust shells and characterize the physical distribution of dust. All the objects were observed 
with the MIDI instrument at the Very large Telescope Interferometer (VLTI) in the 8-13 fim 
wavelength range, typically with a baseline of ~ 47 m. NaStl and WR 106 were also observed 
at 10.7 /xm with a multi-aperture interferometric technique using individual segments of the 
Keck-I telescope mirror to get complementary data from baselines up to the 10 m diameter 
of the Keck mirror. In Sections 2 and 3, we describe the observations, data reduction and 
present the results. Section 4 details the dust models we have constructed and the comparison 
with existing models based on the spectral energy distribution. Section 5 summarizes these 
results and presents the conclusions. 



2. Observations and Data Reduction 
2.1. MIDI 

Most targets were observed over four nights in 2004 July, except for AG Car which was 
allocated one night in 2004 April (Table 1). NaStl was also observed on two nights in 2004 
September. All observations were with the ~ 47 m long UT2-UT3 baseline, except for the 
2004 Septe mber observations of NaStl which used the ~ 89 m UT2-UT4 baseline. The MIDI 



instrument (ILeinert et al.ll2003l ) on the VLTI combines the light from a given pair of 8.2 m 



telescopes (UTs) and produces spectrally dispersed fringes in the 8-13 /im ban d. A detailed 



descri ption of the observing technique and data reduction can be found in iLeinert et al. 



(I2OO4J ). Here we limit ourselves to a brief summary. 



The observing sequence starts with acquiring the target on each individual telescope 
with the 8.7 fim filter (1.4 fim wide) and attaining maximal overlap of the two images, while 
chopping the secondary at ~ 2 Hz over 10". The images are combined to produce the two 
complementary interferograms and dispersed through a prism with a resolution of X/6X ~ 30. 
The delay lines are scanned to find and track the fringes and fringe visibilities are recorded 
while stepping the delay through ~ 10 wavelengths (A =10 fim) in steps of 2 fim with a 
piezo-driven mirror. The scan rate is tuned to the timescale of instrumental and atmospheric 
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fluctuations of delay and a typical measurement records ~ 200 scans in about 5 minutes. 
After the fringe measurement, a photometric measurement is carried out by blocking the light 
from one and then the other telescope, while chopping the secondary mirror. No chopping 
is done during the fringe-tracking scans. The pointing is maintained by the Coude guiding 
with only tip-tilt correction; higher order AO was yet to be integrated into the system. 
A sensitivity better than ~ 1 Jy is achieved for this mode of fringe measurement. Given 
that these observations were during the first regular semester after MIDI commisioning, the 
accuracy that we achieve for the visibility measurements is of the order of 10% and is mostly 
set by photometry, guiding (overlap of the two images) and fringe tracking errors which 
impose limits on the accuracy of calibration. We calculate the errors on the visibility as the 
standard deviation over the multiple calibrations for each measurement and conservatively 
quote the higher of this number or 10% of the visibility amplitude. The above is true for the 
absolute visibility amplitude errors; the relative errors are less and the slope of the visibility 
in the 8 to 13 /im band is accurate to ~3%. 

Standard calibrator stars of known diameter from the MIDI catalogu^ were observed 
immediately before and after the target (Table 1). Given the relatively large systematic 
errors inherent in the visibility determination, we have also used calibrators from other 
programs within a span of 2-5 hours (depending on atmospheric conditions). 

The visibility reduction me thod generally follows the now w ell-established algorithms for 



optical/IR interferometry (e.g., ICoude Du Foresto et al.l (119971 ) for fiber stellar interferome- 



ters: the basic principles apply for non-fiber systems like MIDI) and was carried out using 
custom software provided by the VLTI consortium written in the IDL language. The spectra 
are extracted with a typical binning of 4 pixels, with each pixel being ~ 0.05 fim wide. The 
two complementary outputs of the combiner are differenced to remove the incoherent flux 
including overall variations in the contribution from the two telescopes. Next, the fringe 
power is estimated from the Fourier transform and the visibility amplitud^, normalized by 
the photometric fluxes, is calculated. The visibility is then calibrated using the instrument 
visibility measured on the calibrator stars. 



^ [[pREFjhttp: / /www. eso.org/ projects / vlti/instru/midi / midijcalib.txt 

^The usual estimator used in optical/infrared interferometry is the square of the visibility amplitude 
(V^) since this quantity can be corrected for detector noise bias. However the MIDI measurement error in 
the mode we used ("High-Sens", with no simultaneous photometry) is dominated by calibration error, not 
detector noise. 
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2.2. Keck 



We carried out observations at 10.7 /xm with the Keck-I telescope on two stars in our 
sample, WR 106 and NaStl, in 2005 May (Table 2]. These observatio ns were carried out 
in a "segment-tilting" mode (IMonnier et al.l l2004bl : IWeiner et al.l |2006| ) to achieve close to 
diffraction limited resolution. Imaging at 10 /^m with the full pupil has, in practice, failed 
to reliably yield diffraction limited images. This is in part attributable to the "seeing spike" 
problem, where the visibilities at shorter baselines are difficult to calibrate in the presence of 
changes in the seeing between the target and calibrator. The results shown here demonstrate 
the utility of this new method in alleviating some of these problems. In this technique, custom 
software is used to reconfigure the segmented primary mirror so as to allow selected groups of 
segments to focus on the detector (the Long- Wavelength Spectrograph camera). Each group 
(four in our case) of 6 segments is chosen so as to form a non-redundant sparse-aperture array. 
The 1.8 m hexagonal segments are individually phased up so that each group focuses on a 
separate spot on the detector. Thus in a given frame, 24 out of 36 segments are used in these 
four "Fizeau arrays" (the rest of the segments are pointed away from the detector). This 
arrangement measures 15 baselines ranging from about 2 m to 9 m and 20 closure phases. 
Each frame is of the order of 10 ms and we co-add 9 frames, giving an effective exposure 
of 90 ms. Chopping is done at a 5 Hz rate. Calibrators are observed before and after the 
targets and each target is visited at least twice during a night. The analysis of the data 
closel y follows the apert ure- masking analysis described in detail in numerous publications 



Tuthill et al.ll2000l ). These measurements complement the long-baseline VLTI data by 



(e.g., 

adding in visibilities at shorter baselines up to 9 m. 



3. Results 

3.1. The WC stars 

The two WC type stars in our sample are WR 95 and WR 106. They are both classified 
as late type (WC 9) and are known to be p ersistent dust producers from IR excess detected 



over long-term photometry (IWilliamslll995l ). The photometric data do not have any period- 
icity and neither do the spectra indicate any clear signs of binarity. We obtained two MIDI 
visibility measurements on WR 95 on consecutive nights in July 2004. Though the projected 
baseline and position angle differ slightly between these two measurements, the visibility 
values were remarkably similar (Figure 1, broken lines). Though this might indicate a small 
degree of asymmetry, the accuracy of these measurements are not high enough to justify 
this interpretation. Instead we have chosen the average of the two sets to represent the 
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spectrally dispersed visibility (Figure 1, solid line). As can be seen in Figure 1, we clearly 
resolve the WR 95 dust shell with visibility values ranging from ~ 0.2 to ~ 0.6 over the 8-13 
/im band. In the absence of detailed knowledge of the source geometry we have assumed a 
gaussian profile and the resulting FWHM in milliarcseconds (versus wavelength) is shown 
by the dotted line. WR 95 was not observed with the Keck segment-tilting method since it 
was below the sensitivity limit. 

Figure 2 shows the MIDI results for WR 106. Though we obtained two measurements 
in July 2004 (see Table 1), poor seeing conditions resulted in variations in fringe tracking 
and calibration between the two observations. In the absence of a clear criterion to select 
one data set over the other, we average the two and the resulting larger error bars are shown 
in the figure. However, this target is well-resolved and meaningful conclusions can be drawn 
even from the less-than-optimum data. The visibility variation with wavelength is similar to 
that for WR 95. 

WR 106 was resolved by the Keck single aperture (Table 2). This was surprising in that, 
given the visibility values measured by the VLTI, one would have expected an unresolved 
source at the comparatively short baselines offered by the Keck mirror. This indicates that 
there is an extended component in the dust distribution which was over-resolved and hence 
failed to show up in the VLTI observations. We present these results in Section 4 where we 
discuss some simple models. 

In Table 3, we list the angular sizes (gaussian FWHM) measured with MIDI at 10.5 
/im and the corresponding linear sizes of the objects, assuming distances from the literature 
(see table footnote). The parameters directly constrained by this data set are the size of 
the dust shells and its variation with wavelength. The radius at which the dust grains start 
to conden se is a much - debated and cr ucial topic in the theory of dust-formation around 



WR stars (jZubkdll998l : ICrowtherll2003l ). The major difficulty has been that the equilibrium 
temperatures at the radii predicted by the SED-based models are often above the sublimation 
temperature for carbon dust. Dust formation in the shock interface of a colliding- wind system 
has found some favor in this context. A recent sur nmary of results frona a long-term campaign 



of aperture- masking interferometry at the Keck (iMonnier et al.l 120071 ) provides the near-IR 
sizes for a number of WR stars, including those in our sample. The authors point to a 
tight correlation between the the surface-brightness in the 1.5 - 3.0 fim band and the sizes 
(gaussian FWHM) of the objects in their sample to argue for a common underlying dust- 
production mechanism. Our direct measurement of the mid-IR size-scale for these objects 
is an important addition to this debate. From our gaussian FWHM curves, it is clear that 
the size of the emitting region is more or less invariant in the 8-13 /im band. This could 
indicate that the emission is dominated by an inner edge or ring of material. In fact the sizes 
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for WR 95 at the shorter wavelengths (1.6 and 2.2 fim) reported in iMonnier et al.l (120071 ) 
are quite similar to our values. For WR 106, the mid-IR size is larger, as expected from the 
presence of the extension seen in the Keck measurement. In section 4, we present radiative 
transfer models for the dust shells based on our measurements and derive inner edge radii 
which we compare with values from SED-based models. 



3.2. NaStl 



NaStl has a c heckered history of classif ication. It was initially typed as a late-type 
WN star, WR 122 (Ivan der Hucht et al.lll98l[), but subsequen tly this classification has been 
questioned. The most recent study (ICrowther fc Smithlll999l ) has resulted in this star not 
being considered a Wolf-Rayet any longer. Crowther and Smith find the chemical compo- 
sition of the optical nebulosity sorrounding the star to be similar in many respects to an 
LBV, like rj Carina, and conjecture that the obscured central star could be an early type 
WN, presumably just leaving the LBV phase. The uncertainty in classification is in part 
due to the rarity of these very massive stars, especially one caught in transition. As such 
our measurement of the extent of the close-in dust and its variation with wavelength is of 
significance. The source is unresolved by MIDI in single dish imaging mode and the flux 
measured within the ~300 mas beam is similar to the the value quoted in Smith and Houck 
(2001, Figure 4 in their paper) from a bigger beam, indicating that all the mid-IR emission 
originates from a compact object. Another notable feature is the absence (as in the WCs) 
of any silicate feature in the visibility profile (Figure 3) or the MIDI spectrum (Figure 4). 
Although this does not rule out the presence of silicate dust, it attests to the fact that there 
are no structures w ith silicate emissiori at the MIDI resolution scale. This is in agreement 
with the finding of ICrowther fc SmithI (119991 ) that the nebula is highly enriched in CNO 
processed material. 

NaStl was not resolved with the Keck single aperture. In our experience, the segment- 
tilting method is capable of resolving sizes down to 40-50 mas at the 3-sigma level. At 10.7 
/xm, this is approximately a factor of 6 better than the formal diffraction limit (1.22 A/D) for 
the 10 m aperture. A more rigorous analysis for the resolution limit for aperture-masking 
observations at Keck (the near-infrare d analogue of segrn ent-tilting) estimates it to be 4 
times better than the diffraction limit (IMonnier et al.l 120071 ). We estimate an upper limit of 
50 mas for NaStl from the Keck experiment. 

NaStl was observed with MIDI at both the UT2-UT3 (projected baseline: 44.0 m, 
PA: 42°) and the longer UT2-UT4 (projected baseline: 63.0 m, PA: 78°). Figure 3 shows 
the visibility values and corresponding gaussian sizes for NaStl. The visibility varies with 
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wavelength, decreasing towards the red end of the spectrum. This typically indicates cooler 
material further away from the star and is in contrast to the behavior we see in the WC stars. 
At least part of the dust around LBVs could be swept from the ISM and/or formed in a 
cooler environment in the episodic ejecta from the LBV phase and could lead to a less-defined 
inner edge and more spread in the dust temperature. The visibility on the PA 42° shorter 
baseline is very similar to that on the longer one, and if we intrepret thi s directly as a size, 
may indicate an elongation in this direction. The optical image of NaStl (jCrowther fc Smith 



19991 ) in the Nil A6583 line shows an elliptical nebula with the major axis at a PA of 30°and 
an aspect ratio of 0.6, though the size scale (~ 8 arcsec) is several orders of magnitude larger 
than the milliarcsecond scales we measure. 

However, it must be stressed that the visibility curves from the two different baselines 
are almost identical and the probability of obtaining the same profile for a spatially complex 
object (even if as simple as an ellipse) by pure ch ance is very low. A possib le alternative 



interpretation is that, following the suggestion of ICrowther fc SmithI (119991 ). the mid-IR 
source could be split into two components: a compact, unresolved part, with a correlated 
flux of about 5 Jy from 8 to 13 fim that corresponds to the central source and perhaps some 
hot dust in its immediate surrounding, and a detached rj Car-like nebula, fully resolved by 
the interferometer. The extension of this nebula may be of the same order as that of the 
region contributing the [Nell] 12.8 /xm line which is seen in the single dish spectrum, but not 
in the correlated (i.e., arising in the unresolved core) flux (Figure 4). However, the upper 
limit on the size of this component would be limited to ~ 50-60 mas at most; anything 
bigger would have been resolved by the Keck observations. Given the 10% error on the 
MIDI visibilities, this could still be over-resolved, though marginally, even at the shorter 
44 m baseline. The size of the compact component is possibly less than 5 mas or so, given 
that higher values would start to affect the relative slope of the visibilities between the two 
baselines. Qualitatively, this hypothesis offers a natural explanation of the measurements: 
the central source is unresolved and the extension overresolved whatever the baseline length 
(44 m or 63 m). However, the tight constraints on the extended component weaken the case. 
Clearly, more observations are needed to clarify th e picture. A very h o t sour ce (T~100,000 



K) is required for the fitting the object's SED by ICrowther fc SmithI (119991 ) . and the 5 Jy 
(1.5 X 10^^^ Wm~^yum~^) correlated flux we measure at 10 fan is not consistent with this. 
A constant correlated flux from 8 to 13/im (Figure 4) can hardly be attributed to a source at 
a single temperature (say, ~ 2000 K) either. In order to decompose the SED into multiple 
sources at different temperatures it is important to study this object in the near-IR using 
the AMBER instrument at the VLTI to resolve the K band continuum and the very strong 



Hel 2.058 /zm line that may originate in the close vicinity of the star (jCrowther fc Smith 



19991 . their Figure 6). 
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At a distance of 2 kpc, the maximum linear size from the gaussian model is ~ 40 AU 
and with a two-component model it is ~ 100 AU (for comparison, th e average extension of 
the optical nebula is ~ 13600 AU a t 2 kpc from ICrowther fc Smith! ) . If we take the same 
characteristic expansion velocity as ICrowther fc SmithI (15kms~^) then the dynamical age 
of the dusty nebula at 100 AU from the source is about 30 yr. This shortens to a couple 
of year s if instead we assume dust grains fo rming in the wind at a velocity of ~ 100 -300 
kms~^( 1Zubkdll998l : ICrowther fc Smithlll999l ). Hence the dust that we resolve points to very 
recent or ongoing dust formation close to the star which has implications for its current 
poorly-known evolutionary state. 



3.3. AG Car 



AG Car is one of the best studied LBVs and a prototype of its class. AG Car has 
been extensively observed from UV to radio wav elengths with photomet ric, spectroscopic, 
and polarimetric techniques over several decades (IHumphreys et al.lll989l ). which locates it 
in the H-R diagram very close to t he Humphreys-Davidson limit, an empirical luminosity 
boundary for evolved massive stars (IHumphreys fc Davidson! Il979l ). 



AG Car is surrounded by a bright extended bipolar nebula, 3 0" x 40" in size fjNota et al. 



19921 ). perhaps from a massive erruptive event 10^ years ago flLamers et al.ll200lh . The 



kinematics of the bipolar nebula is consistent with a two-lobe polar ejection view ed nearly 
from the equator (INota et al.lll992l ). Spectroscopic analysis by iGroh et al.l (120061 ) supports 
the bi-polar geometry and shows the 50-100 R0 supergiant central star to be rotating close 
to its break-up velocity. The high rotational velocity after the huge LBV outburst may be 
understood if the original explosion expelled the mass and momentum preferentially toward 
the poles. The wind at smaller latitude may not have formed or may have launched material 
with velocities lower than the es cape velocity whic h fell back onto the star. The nebula 



contains a large amount of dust (IVoors et al.l I2OOOI . ~0.2 Mq), which implies a huge mass 



loss rate during the eruption (at least 8-15 M, 



Our measu r emen ts are set in con text by two s tudie s in particular; IVoors et al.l (I2OOOI ) 
and iNota et aD J2002h . The images in IVoors et aU koO(h obtained with the 10 /im TIMMI 
camera on the 3.6 m NTT at La Silla reveal a central peak and a detached nebula with a 
radius of ~ 10". Their SED analysis reveals a complex situation: the dust is dominantly 
oxygen-rich but Polycyclic Aromatic Hydrocarbon (PAH) features are observed at 3.3, 7.7, 
8.6 and 11.3 /im revealing some amount of carbon dust. Moreover, a population of very 
large grains (> lO/im) is necessary to explain the flux levels at long wavelengths whereas a 
population of small warm grains is invoked to explain the flux between 5 and 20/im. 
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There is evidence for CO J = 1^0 and 2^1 emission associated with AG Car (INota et al. 



2OO2I ) and the authors conclude that the unresolved CO emission must arise close to the star 
in a high density warm region and not in the detached nebulae resolved in optical emission- 
line images. They postulate that a warm and dense equatorial disk exists close to the star. 

MIDI observations of AG Car were carried out with the UT2-UT3 baseline on two 
nights, with projected baselines of 42.0 m and 35.6 m and PA of 41° and 69° respectively. 
The visibilites and corresponding gaussian sizes are shown in Figure 5 and the spectrum in 
Figure 6. AG Car is marginally resolved. The measu red size (~ 11 mas) is similar for both 



observations and fairly uniform with wavelength. The lNota et al.l (|2002| ) optical images show 
strongly enhanced Ha brightness at PA~35-225°, and the morphology is similar at 12/im 
(IVoors et al.ll2000l ). This direction is also the waist of the nebula, the largest extension being 
seen perpendicular to it. The MIDI baselines are therefore oriented close to an optimal 
direction to detect equatorial material and we are very likely resolving the equatorial disk, 
though we cannot presently distinguish between a disk and shell structure. This is the first 
direct evidence for dust in a compact structure. The size we detect (11 mas) is ~ 3 orders 
of magnitude s malle r than the upper limit (3") to the size of the (unresolved) CO outflow 



m iNota et al.l (120021 ). The corresponding linear size is of the order of only 100 AU at a 
dist ance of 6 kpc; ver y close to the central star. Nevertheless, this detection adds weight to 
the iNota et al.l (120021 ) claim that the CO emission arises from molecular gas close to the star 
and is shielded by a disk s tructure. The pre sence of dust has implications also for the origin 
of this compact structure. iNota et al.l (120021 ) f avor a scenario where the the star goes through 
a Red Super Giant (RSG) or RSG-like phase (ISmith et al.lll998l : iLamers et al.ll200ll ) during 
which the m olecular outflow is formed. Dust format ion is also commonly observed during the 
RSG stage (IMassey et al.ll2005l : iDanchi et al.lll994j ) and fits in with this hypothesis, though 
the ~ 100 AU scale size indicates very recent or current dust formation as well. 

There are few published spectra of AG Car in the mid-IR and the MIDI spectrum is a 
valuable addition. For comparison we also show in Figure 6 the spectrum from the Short 
Wavele ngth Spectromete r (SWS) on the Infrared Space Observatory (ISO) in the 8-13 /xm 
range (ISloan et al.l l2003l ). scaled to roughly match the flux level of the MIDI spectrum. 
The ISO spectrum, with a much bigger aperture, exhibits a slightly different slope from the 
MIDI measurement, but the striking difference is that the MIDI spectrum does not show the 
11.3/im PAH fe ature, imply ing that the PAH originates from further out in the disk/shell as 



conjectured by IVoors et al. 



- 11 - 



The WC Models 



The dust around late type WR stars has been closely studied over the last two decades 
or so and presents a challenging problem pertaining to its formation, geometry and evolution. 
There have been quite a few attempts to model the dust in some detail. WvdHT model these 
as spherical shells based on extensive IR spectroscopy. IZubkd (Il998l ) provides dust models 
with emphasis on the physics of grain forma ti on an d growth taking into account the grain 
dynamics in the stellar win d. iHarries et al.l (12004r) take advantage of the high resolution 



near-IR images of WR 104 (ITuthill et al.lll999l ). which show a "pinwheel" of dust formed 
in a binary wind interface, to construct a 3D radiative transfer model for this object. As 
mentioned earlier, a large percentage of the known dust-forming late type WC stars are 
known to be binaries, WR 104 being a prime example. The forma t ion of dust in these stars 
is usually interpreted using the wind-wind collision model of lUsovl Jl99lh. None of the stars 
in our sample show definite evidence of being binaries (jWilliams &: van der Huchtl |2000| ) , 
though this remains a strong possibility. 

It has been difficult to explain dust formation in the hostile environment of the wind 
from single WR stars. SED-based mo dels typically have dust forming at distances from 300 
up to a 1000 times the stellar radius (jZubkdll998l ). The consensus on the material of the 
dust, from the SED signature, is amorphous carbon (WvdHT). However, at t hese distance s 
the electron temperatures in the winds of WRs are likely as high as 10^ K (jZubkd Il998l ). 
implying that carbon is almost fully ionised. 

We present here radiative transfer models for the WC stars in our sample with an aim 
to use the interferometric measurements to answer the crucial question as formation radius 
of the inner edges of these dust shells. To this end, we use a 2D Monte Carlo radiative 
transfer code, MCTRANSF (INiccolini et al.ll2003l ). to construct spherically symmetric dust 
distributions consistent with the measured visibilities. 

The parameters of the models are the star diameter and effective temperature, the 
internal and external radius of the envelope, its density profile and optical depth, the minimal 
and maximal radius and size distribution of the grains and the optical indices of the dust. 
We start with the basic parameters for the dust shells around WR 95 and WR 106 mostly 
taken from the comprehensive SED-based models of WvdHT (see Table 4). By propagating 
energy packets in the sampled envelope (circular grid), taking into account the absorption, 
emission and scattering characteristics of the dust, MCTRANSF computes the temperature 
of the envelope, the SED of the flux received by the observer directly from the star, as well 
as the flux emitted and scatt ered by the envelope. The code has been described in detail 



m 



Niccolini et al.l (120031 ) (see iPascucci et al.l . |2004| . for benc hmarking tests). An exam ple 



application in modeling a circumstellar dust shell is found in lWoitke fc Niccolini! (120051 ). 
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We then explore the parameter space around this model by primarily varying the lo- 
cation of the inner edge of the disk to see what values are most consistent with both our 
measured visibilities and spectra. We also tried limited variations of the optical depth, grain 
size and distribution to quantify their effect on the visibility profile and spectra. The ra- 
dius of the central star and the composition of the dust were not changed in these trials. 
MCTRANSF does not explicitly incorporate the sublimation temperature of the grain as an 
input. The dust temperature in the nucleation z one from ex i sting observational and model- 
ing constraints ranges from 1000 - 1500 K or so (jZubkdll998l : iHarries et al.l 12004 ) and in the 
process of iteration we discard those models which produce untenably high temperatures at 
the inner edge of the shell. We stress that the modeling is not a "fitting" exercise. That does 
not seem justified, given the limited visibility data and the fact that the actual geometry of 
the dust is most likely more complex than simple spherical shell models. Our goal here is to 
use the new measurements to constrain the vitally important distance to the dust-forming 
zone from the star. 



4.1. WR 106 



Extensive studies of the SED of this well known dust- producing WC9 star exist i n the 
literature and long-term IR photometry is summarised in Ivan der Hucht et al.l (120011 ) . In- 
frared (J-Q bands) measurera e nts ar e presented in WvdHT. More recent mid-IR ISO spectra 
are found in ISmith fc Houckl (120011 ). Given the difficulty in standardising mid-IR photom- 
etry from various sources and the fact that dusty WC stars can have differing degrees of 
variability, we rely mostly on t he 8-13 /im spectra o btained with the MIDI instrument simul- 
taneously with the visibilities. iLeinert et al.l (120041 ) have confirmed the reliability of spectra 
from MIDI. We used standard reduction techniques as they describe to extract the spectra; 
most of our interferometric calibrators were also chosen to be IR spectroscopic standards □ 
and as such provide ideal templates to calibrate the spectra. De-reddening WR spectra has 
been non-trivial, primarily because these stars are grouped towards the Galactic center with 
large columns of intervening interstellar dust in addition to the circumstellar dust. The wide 
interstellar silicate feature centered around 9.7 /im comphcates the issue ( ISmith fc Houck 
200ll ). To provide a uniform comparison, we have de-reddened the 8-13 /xm band with 
values of Aa/Av from WvdHT and the van der Hulst extinction curve 15 for the shorter 
wavelengths. 



For WR 106, the model visibilities from our "best" model are shown in Figure 7, along 



^TlIR.EF]http://ssc.spitzer.caltech.cdu/irs/calib/tcnipl/cohen_models/ 
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with the Keck measurements as well as the MIDI visibility at 10.5 fim. 

The MIDI visibilities across the entire band (8-13 yum) are compared with the model 
in Figure 8. The closure phases measured by the Keck segment-tilting experiment did not 
show any significant deviation from zero, signifying a lack of asymmetry at this resolution 
and sensitivity. We therefore plot the Keck visibilities treating the baselines as scalars (i.e., 
ignoring the orientation). Figure 9 shows the model SED along with our measurements 
and those from the literaturcl spanning wavelengths from the optical to the far-infrared for 
comparison. The temperature profile of the dust in our model is shown in Figure 10. The 
parameters for the model are as below: 



Central Star 



For the central star, we have adopted a blackbody spectrum at the effective temperature 
listed below. The temperature of the WR star has been assigned a range of values in the 
literature. Most dust-formation models adopt a value between 18000 and 25000 K (e.g., 
WvdHT). These values have bee n criticized as being under-estimates from recent studies of 
the WR spectra (ICrowtherll2003l ). However, we find in our models that higher temperatures 
(we tested up to 40000 K) do not fit the 8-13 /xm spectra and also lead to dust temperatures 
above the 1500 K or so sublimation temperature usually adopted for amorphous carbon 
grains. This indeed has been the reason that previous modeling attempts favour the lower 
values as well. Although more sophisticated SEDs have been used for the stellar photo- 
sphere (IHarries et al.ll2004j . WdvHT), a simple blackbody is likely adequate for our purpose. 
The radius of the st ar is is assumed to be I4.6R0 (WdvHT) at a distance of 2.09 Kpc 
Jvan der HuchtlboOlh . 



* Temperature of central star (Tg//) = 19000 K 

* Angular radius of central star (R,,) = 3.23x10^^ arcsecond 



Dust Shell physical properties 



Composition: Amorphous carbon (optical constants from lRouleau fc Martini . 1991 



''From catalogs available at the VizieR online database (jOclisenbein et al.ll200Cll ). Optical fluxes: Homo- 
geneous Means in the UBV System (Mermilliod, 1991); Hipparcos Input Catalogue, Version 2 (Turon et al., 
1993); The USNO-Bl.O Catalog (Monet et al., 2003); IR fluxes: Catalog of Infrared Observations, Edition 5 
(Gezari et al., 1999); Third release of DENIS (DENIS consortium, 2005). 
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* Grain size: Ranges from 0.5 to 0.6 fim 

* Optical depth (r) = 0.015 at 10 fim 



Dust Shell geometry 

* Inner radius of dust shell (Rmt) = 280 R* (multiples of stellar radius) 

* Outer radius of dust shell = 8400 R=^ 

* Density profile: Proportional to r~^, except in the region between 1.0 Hint and 14.0 
Rint where we introduce an over-dense region with density proportional to 7xr~^. 



4.2. Discussion 



Most of the starting parameters of the dust shell model were based on WvdHT. They 
establish the composition to be amorphous C at least in the case of the episodic dust-producer 
WR 140. The lack of silicate features in the mid-IR region also points to amorphous C. The 
grain size and distribution for the shells are not well known. The Zubkol ( 1998) grain- 
forma tion model predicts a fairly small maximum grain size(~ 100 - 200A). Harries et all 
( I2OO4J ) model the inner part of the WR 104 pinwheel nebula with small (~ lOOA) grains, 



mainly to match the 2-10 fim flux. However, t here is evidence in the SEP of WR neb ulae 
pointing towards l arge g rains of ~ 0.1 to 1 fxm Jchiar fc TielenslboOll : IVeen et allEggsl ) and 
Marchenko et al.l (120021 ) find a characteristic size of 0.5 fim for the extended emission in WR 
112. The model presented here also favors large grains with a quasi-uniform grai n size ranging 



from 0.5 to 0.6 /im, with a standard interstellar distribution scaling as a ^'^ (jMathis et al. 



19771 ). A model with sizes ranging from 0.03 to 2.0 /imalso shows reasonable agreement with 
the data for similar values of the other parameters. Despite extensive trials we were unable 
to find a distribution of smaller grain sizes which agreed with the data, in particular the fiux 
at the shorter infrared wavelengths. 



The optical depth for these dust shells in the mid IR is fairly low (IZubkdll998l ). We use 
here a value of 0.015. Our primary goal of matching the measured visibilities is not strongly 
affected by the optical depth and it is the spectra that constrain this parameter. 

The important constraint imposed by the visibility profile is on the range of inner edge 
radii of the dust shell. The "best" model in Figure 7 has a value of 280 R* for the inner 
edge radius, which is the same as in the WvdHT model. We show in Figure 8 the model 
visibility profile for values ranging from 230 R* to 370 R*. Even for this narrow range, the 
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visibility profiles are clearly different, indicating the fairly strong constraints imposed by the 
measurements. 

The over-density (denser than the r~^ profile) that we introduce in the shell is to ac- 
comodate the visibilities measured with the smaller baselines at the Keck, which are lower 
than expected for a r~^ profile consistent with the longer basehne MIDI visibihties. The 
over-dense region in the dust shell extends from the inner edge (Rmt) up to 14 Hint in the 
model, i.e., size scales from 9 to 125 mas. This enhancement at the spatial scales probed 
by the segment-tilting experiment (resolution ~ 50 mas, see Section 3.2) decreases the vis- 
ibility at the low spatial frequencies. This of course is not a unique solution, but offers an 
intuitive physical model. The over-density could indicate an episode of more vigorous dust 
production. It could also be an artifact of applying an overly simplistic spherically symmet- 
ric model to a possibly complex shell geometry as in, for example, a "pin-wheel" structure. 
In the absence of more complete visibility sampling, we limit ourselves to simply pointing 
out the existence of an extended structure in the dust around WR 106. 



4.3. WR 95 

WR 95 is also a late-type dust-forming WC, similar to WR 106. As in the case of WR 
106, there is no definite evidence so far, either spectroscopic or photometric, of binarity. WR 
95 falls below the sensitivity hmit of the Keck segment-tilting experiment and so we base 
our models solely on the MIDI visibilities and spectra. Our best model visibilities along with 
the measured values in the 8-13 /xm MIDI band are shown in Figure 11. Figure 12 shows 
the measured 8-13 fim SED, the near-IR and far-IR photometry from the literature, and the 
model SED profile. The dust temperature profile is as in Figure 13. The parameters of the 
model are as below: 



Central Star 

* Effective temperature (Te//): 25000 K 

* Angular radius of star (R*): 3.2x10"^ arcsecond 

Dust shell physical properties 

* Composition: Amorphous carbon 
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* Grain size: Ranges from 0.5 to 0.6 fim 

* Optical depth (r) = 0.005 at 10 /xm 

Dust Shell geometry 

* Inner radius of dust shell (Rmt) = 480 R,, (multiples of stellar radius) 

* Outer radius of dust shell = 1230 R* 

* Density profile: Proportional to r^^ 

As in the case of WR 106, we base our initial model parameters mostly on WvdHT. To match 
the near IR photometry and our visibility measurements we had to adopt higher temperature 
of 25000 K for the WR star. The optical depth for our best model, 0.005 at 10 /im, is lower 
than that for WR 106. Since we do not have shorter baseline Keck observations in this case 
to detect extended structure, a uniform r~^ density profile suffices. The nominal inner radius 
of the shell is 480 R* (Figure 11 shows how the model visibilities vary for slightly different 
values of this parameter), compared to 410 R* derived by WvdHT. As in the case of WR 
106, a dust grain distribution with sizes ranging from 0.3 to 2.0 fim also yields a fair fit to 
the measurements. 



5. Conclusions 

High-resolution measurements of dust around massive stars are sparse in the literature. 
Coupled with the rarity of the objects themselves, it is not surprising that so much remains 
unclear about the mechanisms of dust-formation and the nature of the dust itself. Only a 
handful of WR stars known to produce dust have been imaged with any degree of detail. 
Three of these, WR 104, WR 98a and WR 140 are binaries, producing dust at the wind- 
wind interface. The fourth, WR 112 reveals a multiple- shell dust structure, though at 
much l arger (arcsecond) scales than discussed here for our targets, suggestive also of being a 



binary (IMarchenko et al.ll2002l ). Hence the measurements presented here, albeit with limited 
baselines, of similar objects with no previous evidence of binarity, are an important step in 
probing the nature of these dust shells. 

We observed two WC9 stars (WR 106 and WR 95), an LBV (AG Car) and an object 
which may be in a stage of transition from an LBV to a WR phase, NaStl. NaStl was well 
resolved and the visibility measurements suggest increasing size with wavelength, a behavior 
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unlike the WC 9 stars in our sample. We favor a two-component model for the dust with a 
compact, unresolved core and an extended component. Previous studies of the SED of this 
object have also indicated this structure. For AG Car as well, our measurements resolve a 
compact mid-IR source which is likely associated with a known molecular gas component 
and is the first direct evidence for dust near the central star. For both these objects, this 
study has revealed the presence of close-in dust that has important implications for the state 
of evolution of the central star. 

For the WC stars we construct spherically symmetric radiative transfer models for the 
dust shells. These models indicate the inner edge of these shells are at a few hundred 
stellar radii (tens of AU) from the star. Gaussian sizes (FWHM) directly estimated from 
the visibility measurement provide model-independent upper limits which are consistent with 
these values. At these distances, the dust temperature is fairly close to sublimation values for 
carbonaceous dust. In addition, our results support previous indications of fairly large grains 
(~ 1 /xm) despite these conditions. These new results re-open an existing controversy: How 
does dust form in the hostile environment of a WR stellar wind at such high temperatures? 
The proximity of the inner edge to the star also forces us to adopt values close to 20,000 
K for the central star to keep the dust below the sublimation point. However, the spectra 
of late-type WC stars point to higher (40,000 to 70,000 K) temperatures, rendering the 
situation even worse. Our results certainly indicate that dust formation is likely occuring 
under fairly extreme conditions. This may suggest a coUiding-wind binary model but direct 
evidence is still lacking for the objects in our sample. 

The dust around WR 106 seems to deviate from an r"^ profile at the shorter baseline 
Keck measurements. Our simple models suggest an over-dense region fairly close to the inner 
edge. This could be also be a first indication of of binarity with multiple shells from episodic 
dust production or a more complex geometry as in a pinwheel. But given our limited baseline 
coverage, a caveat on modeling artifacts is warranted. 

We have not attempted to estimate the dust mass from the models. Our measurements 
are insensitive to the colder dust further from the star and it is unlikely that we can improve 
on the existing dust mass estimates for these stars from spectral profiles. 

In the long-term, further progress would involve multiple baselines to image the dust 
structure. Observations of the outer dust at longer wavelengths are also required to accu- 
rately estimate the dust mass and geometry. In the short-term, we aim to look for asymme- 
tries in the dust distribution by using the UT3-UT4 baseline at the VLTI, orthogonal to the 
one used here. 

This work was performed in part under contract with the Jet Propulsion Laboratory 
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the University of California and NASA. The Observatory was made possible by the gen- 
erous financial support of the W.M. Keck Foundation. The authors wish to recognize and 
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Table 1: MIDI Observations 



Star 


V 
mag 


IRAS U/im 

Jy 


lype 


Date Observed 


Calibrators 


WR 95 


13.4 


4.56 


WC 9 


07-30-2004 


HD 177716, HD 129456, HD 169916 










07-31-2004 


HD 169916, HD 165135, HD 192947 


WR 106 


11.93 


11.67 


WC 9 


07-09-2004 


HD 165135 


NaStl 


15.4 


14.5 


WN(?) 


07-30-2004 


Same as WR 95 above 










09-29-2004^^ 


HD 168723 










09-30-2004^ 


HD 168723 


AG Car 


7.09 


12.42 


LBV 


04-10-2004 


HD 81797, HD 107446, HD 129456, HD 139997 



"Baseline UT2-UT4 
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Table 2: Keck Observations 



Star 


Date of Observation 


gaussian FWHM 






mas (@ 10.7 /im) 


WR 106 


06-26-2005 


56.0±6.9'' 


NaStl 


06-26-2005 


50.0 



Dominant error arises from calibration 

Object was unresolved. Value is an estimated upper limit 
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Table 3: Measured Sizes (MIDI) 



Star 


Distance 


gaussian FWHM 


Linear Size 


References 




Kpc 


mas (@ 10.5 fim) 


AU 




WR 95 


2.09 


28.4±1.5 


59.4 


1 


WR 106 


2.3 


45.0±5.0 


103.5 


1 


NaStl 


2.0 


20.1, 14.8 (±2.0)^ 


40.2, 29.6 


2 


AG Car 


6.0 


12.5±3.0 


75 


3 



"The two values are from two different baselines. Interpreting sizes from visibility may be inaccurate for this 
object; see text 

References.-fl ivan der Huchtl |200l[) : (2 lCrowther fc SmithI |l999l ): fs lHumphrevs et all |l989l ) 
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Table 4: Starting Parameters for the WC Models 



Parameter 



Value^ 



Teff 

R* 

Dust shell Rint 

Dust shell thickness 

Dust composition 
Density Profile 
Optical depth @ 10 /^m^ 
Grain sizes*' 



1900 K 

14.6 Rq 
410 R* (WR 95) 
280 R* (WR 106) 

3 Rint (WR 95) 
30 Rint (WR 106) 
Amorphous C 
r-2 
0.01 
0.1 to 1 iim 



Unless stated otherwise, all values are from WvdHT 
"Single values are common for both stars 
''See text, Section 4.2, for references 
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Fig. 1. — Left: Visibility profile with wavelength for WR 95. The solid line is the mean of 
the two measurements. Right: The corresponding (mean visibility) gaussian FWHM. 
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Fig. 2. — Left: Visibility profile with wavelength for WR 106 Right: The corresponding 
gaussian FWHM. 



-28- 




Fig. 3. — Left: Visibility profile with wavelength for NaStl (only one set of error bars shown 
for clarity). Right: The corresponding gaussian FWHM. Transforming the visibihty directly 
to a gaussian size is, however, probably inaccurate for this object (see text). 
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Fig. 4. — MIDI spectrum and correlated flux for NaStl. The Ne[II] 12.8 //m feature is not 
seen in the correlated flux and is likely fully resolved by the interferometer as is a large 
fraction (~ 1/2) of the total continuum flux. 
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Fig. 5. — Left: Visibility profile with wavelength for AG Car. Right: Corresponding gaussian 
FWHM. The error bars are similar for both measurements. For clarity only one set is shown. 
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Fig. 6. — MIDI spectrum for AG Car, with the ISO spectrum (scaled overall to match the 
MIDI flux level) for comparison. The PAH feature at 11.3 //m is not seen in the MIDI 
spectrum. 
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Fig. 7. — Visibility versus baseline for WR 106. The solid line is the model visibihty, the 
filled circles are the Keck visibilities (upper left) and the diamond, MIDI at 10.5//m. 
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Fig. 8. — Visibility versus wavelength for WR 106. The line (with error bars) is the MIDI 
measurement. Models with varying inner radii, Rjnt (in units of R*), for the shell are shown 
for comparison. 
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Fig. 9. — SED of WR 106. The solid line is our best model, the grey circles are the 
measurements from MIDI, the squares are measurements from the literature (see text for 
sources) . 
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Fig. 10.~ 



Dust temperature profile for tlie WR 106 best model. 
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Fig. 11. — Visibility versus wavelength for WR 95. The line (with error bars) is the MIDI 
measurement, the filled circles are the visibility points from the best model. For comparison 
we show visibilities for models with different inner radii (Rmt) for the shell. 
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10° 10^ 
Wavelength (/um) 

Fig. 12. — SED for WR 95: The solid line is our best model, the grey circles the measurements 
from MIDI (de-reddened as in text), the diamonds are measurements from the literature (see 
text for sources) 
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13. — Dust temperature profile of our WR 95 model. 



